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Abstract 

The charge dynamics of spin-liquid states described by U{\) gauge theory coupling to fermionic 
spinons is discussed in this paper. We find that the gapless spinons give rise to a power law 
optical conductivity inside the charge gap. The theory is applied to explain the unusual optical 
conductivity observed recently in the organic compound k — {BEDT — TTF)2Cu2{CN)^. We 
also propose an optical experiment to search for the in-gap excitations in the Kagome spin liquid 
insulator. 
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Recent work has shown that the organic compound k — (BEDT — TT F)2C%i2{C N)^£- 
and the spin-1/2 Kagome system ZnCu^OH^Cl^^ hold great promise as the first two 
examples of spin liquid states realized in dimensions great than one^^. Spin liquid states 
are Mott insulators with an odd number of spin | per unit cell which shows no long-range 
magnetic order. They are proposed to exist in systems either in the vicinity of the Mott 
transition 7 or with frustrated lattice structures. In both cases the system can be modelled 
by an appropriate Hubbard model with on-site repulsion U and hopping integral t at half 
filling. For large enough U compared with t, charge excitations are gapped and the system 
is a spin liquid if long-range order is absent in the spin sector. 

The two recently discovered systems are believed to be two-dimensional spin liquids. In 
the case of k — (BEDT — TTF)2Cu2{CN)z the system is described by a Hubbard model on 
a triangular lattice. Since the system can be driven metallic (indeed superconducting) under 
pressure, it is believed that U is not very large compared with t and the insulator is near the 
Mott transition^. In this case the charge excitations acquire a gap, and it is proposed that 
spin-charge separation occurs and spin | excitations (spinons) form a Fermi surface 7 -^ In 
the case of ZnCu^iOH)^^, it is believed that U >> t and the spin dynamics is described 
by the antiferromagnetic Heisenberg model. The frustrated Kagome lattice gives rise to a 
spin liquid state with Dirac fermions excitation spectrum 9 . A unique feature of these spin 
liquid states is that the spin excitations are always coupled to internal U(l) gauge fields 
representing spin-chirality fluctuations 7 - 1 ^ in the spin systems. 

It is often thought the Mott insulator should be fully gaped in its optical (charge) re- 
sponses. Furthermore, the spinons are considered to be neutral and do not absorb elec- 
tromagnetic radiation. Here we point out that due to coupling with the gauge field, the 
spinons do contribute to optical conductivity, yielding a power law absorption at low fre- 
quencies. This may explain some puzzling experimental observations recently reported in 
the organics 10 . 

The dynamics of the spin liquid states can be studied in a slave-rotor representation 
of Hubbard models 7 - with the appropriate lattice structures. In this presentation the elec- 
tron operator is represented as c(c + ) ia = /(/ + )i -e~( + ) j6 ', where f(f + )i a is the spin anni- 
hilation(creation) operator and e~^ + ' lB lowers (raises) the charge by one. The total charge 
operator pi = J2 a furf™ — 1 is the conjugate variable to 6 in this representation. 

After making a mean-field approximation, the low energy effective action of the system 
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can be written in terms of 9 and f(f + ) fields separately, L 1 ^ = L c + L]^, where L c 
representing the charge dynamics and L s represents the spin dynamics of the system. L c is 
described by the strong coupling phase of a quantum x — y modet^, 

i <i,j> 

coupling to internal gauge fields (a , a), where t e ff ~ at with a < 1 being a numerical 
factor determined self-consistently from the mean-field equation, (Aq, A) represents the real 
electromagnetic field coupling to the system and 

l<P = £ " ia » ~ - - 5 ) 2 ^) ( lb ) 

in the case of k — (BEDT — TTF) 2 Cu2{CN) 3 which is believed to possess a spinon Fermi 
surface, [if is the chemical potential, m" 1 is expected to be of order of the exchange 
J ~ t 2 /U. In the case of ZnCu^{OH)^Cl2 where spinons have a Dirac fermion spectrum, 

L i 2) = £ _ l "( a A» + AJ^+a + - + A^T^-c) , (lc) 

//a 

where /i = 0,1,2 and r M are Pauli matrices. The two-component Dirac spinor fields ip± a 
describe two inequivalent Dirac nodes in the spinon spectrum 9 . 

Effects of disorder and phonons can also be included in the actions. Their contributions 
can be included by an adding a term 

to Ls , where V(q) is a disordered potential and b(b + ) 9 are phonon annihilation (creation) 
operators with momentum q. M(q) is the electron-phonon coupling and is the phonon 
dispersion. A corresponding term can also be added to L s for Dirac fermions. 

The thermodynamic and magnetic properties of the above systems have been studied in 
several previous papers^"^. We shall concentrate on the charge dynamics of these spin 
liquid states here. We assume a Mott insulator state with no broken symmetry and with 
isotropy in space. The current response function is given by the conductivity, which can be 
decomposed into longitudinal and transverse parts cry and a±. For a U(l) spin liquid, the 
Ioffe-Larkin composition rule^ relates the physical a to the response function of the spin 
and charge components, 
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<r±(q, w) = (<r B l(q, u) + a c l(q, u)) 1 (2) 

and similarly for cry. Here cr s and cr c are given by the proper response functions of the spin 
and charge (represented by 9) fields appearing in the action L s and L c , respectively. The 
proper response functions represent sum of all diagrams which cannot be separated into 
two parts by cutting one interaction line associated with either the real or internal gauge 
field, and represents the current response of the charges and spinons to the potential a + A 
and a,— ^ respectively. Notice that both the phonon and impurity contributions can be 
included in the definition of the proper response functions. The origin of the Ioffe-Larkin 
rule is that an external A field induces a nonzero a field which is needed to enforce the 
constraint j Cfl + j Sfl = 0.— Thus even though the A field couples only to the 9 field, the 
induced a field indirectly couples to the gapless spinons. 

We parametrize the longitudinal response of the charge field by a dielectric constant e c 
and ignore the analytic correction in q 2 ,u) 2 for small q and uo. Then 

e c = 1 + -■ (3) 

u 

We expect e c — 1 to decrease with increasing charge gap. Furthermore, for small q there 
is no distinction between longitudinal and transverse response in an insulator. Using ([3]) for 
both, we find using Eq. (T5]), 

<T||(JL) (q, w) = "Y • (4) 

w + 4 ^i^rJ <Mi(-q(9>w) 

We should point out that the replacement of the charge response by a dielectric constant 
is not as innocent as it appears. This step should be considered in the spirit of random 
phase approximation and justified using a expansion. The concern is the existence of 
Feynman diagrams involving multiple gauge field lines going across. In the language of 
proper response function, these become part of the charge vertex which couples to the 
external gauge field. Since the gauge field carries gapless excitations, the approximation of 
this vertex by a dielectric constant is not strictly correct except as leading order in -j—-. 

Now we consider the optical conductivity given by cr±(q = 0,u). In this limit there is 
no distinction between longitudinal and transverse and we can drop the _L subscript. The 
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spinon conductivity is expected to be metallic-like. We can safely assume Re[a s (0, u)] » u 
and Im[a s ] << Re[a s ] for small u, and we obtain from Eq. (jlj) 

Re[a(u)] = oo 2 ( £ -^-) " - * Vy (5) 
V 4tt J Re[a s (uj)] 

Note that Rea(u) = for u = as expected for an insulator, but we find contribution inside 
the gap for small uj. First we consider the case when disorder scattering of the spin is weak. 
Then a s (uj) = ne 2 r(u!,T)/m s . The dominant contribution to r _1 is inelastic scattering due 
to the gauge field,— which is given by - ~ [max{fko 1 fceT)] 4//3 . For fko > h/r , kT where To 
is the elastic scattering time, we find 

i^ H ]=u^(^) 2 ^ (6) 

in qualitative agreement with what is observed experimentally in k — (BEDT — 
TT F)2Cu2{C N)^-. Our theory also predicts that Re[a(u>)} crossover to ~ uj 2 at small 
enough u. The above results are strongly modified if localization effect is important and 
ct(uj) vanishes at u — > faster than u. In this case 

o(uj) ~ G s (uj) 

will show similar behavior as observed in usual strongly disordered metals. 

The above analysis can be generalized straightforwardly to the Kagome system 
ZnCu3(OH)§Cl2 which is believed to be a spin liquid with Dirac fermion excitation spec- 
trum u = cq. The only difference is that the "proper" response functions should be replaced 
by the corresponding functions for Dirac fermions. In this case a s (q,u) has the universal 
form 



,,2 (t3.„1 



1+0 



: (7) 

8 too 

where (3 = for non-interacting Dirac fermions and is an unknown exponent in the presence 
of gauge field interaction. Putting <j s (0,u) into Eq. (j4j), we predict Re[a(u>)] oc u 2 " 13 for 
j3 < 1 and Re[a(u)] oc uj 13 for j3 > 1 and we see that the optical conductivity probes directly 
the unknown exponent (3. Since the Kagome system is deep in the Mott insulator regime, 
the observation of power law conductivity inside the Mott gap strong than u A (see Eq. (jl4cj) 
below) will be strong evidence for the existence of gapless spinons. 

We shall now study the general dielectric response e(q,u) of k — {BEDT — 
TTF) 2 Cu2(CN) 3 , which is believed to possess a spinon Fermi surface in more detail. We 



shall assume that the residual interactions are weak enough so that the spinons are in a Fermi 
liquid state. The proper density-density response function of the system is Xd(q, which 
represents the sum of all polarization diagrams which cannot be separated into two parts by 
cutting one Coulomb interaction line associated with the real electromagnetic fielcr^^ The 
dielectric function of the spin liquid is therefore 

e(q,u) = 1 - v e (q)xd(q,uj) 

where v e (q) = 47re 2 /q 2 is the real Coulomb interaction. We assume here that the (3D) 
system is a sum of layers of spin liquid here. 

Charge conservation gives Xd — {q 2 / uj2 )Xh where X\\ is the longitudinal current-current 
response function which is in turn given by a\\ = e 2 x\\/i^- Combining these relations we 
obtain the usual formula 



e(q, uj) = 1 + Ania\\(q, uj)/uj (8) 

where a\\ is related to a s \\ by Eq. (jlj). In the absence of scattering, we expect the density- 
density response function to be 

dn i'juj 

Xds = -T- + — (9) 
an VFq 

where 7 is the quasiparticle density of states at the Fermi level and vf is the Fermi velocity. 
Eq. is valid in Fermi liquid theory and has been shown to remain applicable for small 
q, w when gauge fluctuations are treated to two loop order.— Particle conservation again 
allows us to write a s \\ = iuxds/q 2 - Combining these results we find 
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(10) 



Using Eq. (JHJ), we obtain at small q 



e(q,u)=e c + ^ ^_ } . (11) 

d^i Vpq 

The static dielectric constant is given by the charge part e c and the full dielectric function 
is in principle measurable by electron diffraction. 
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FIG. 1: (la) Leading order Feynman diagram representing coupling between spin and den- 
sity/current fluctuations. Solid lines represent electron Green's functions. There is another dia- 
gram where electron lines reverse in direction, (lb) Corresponding Feynman diagram representing 
correction to proper density-density response function. 

Phonons have small effects on the above results. It only modifies the interaction parame- 
ter 7 and renormalizes the compressibility dn/dfi. The effect of disorder can be included by 
modifying Xdsiq, into a diffusive form ^7 of^+iuj ^ localization effect is not importantr^. 
In this case, we obtain 



where D is the spinon diffusion constant and cr Si || = e 2 ^D. For q = 0, Eq. (fl2|) is consistent 
with the AC conductivity given by Eq. ([5]) as expected. 

It should be emphasized that the coupling of density and current responses to spin exci- 
tations exists rather generally in insulators and does not rely on existence of a spin-liquid 
state. Assuming that the electronic properties of the insulator is described by a Lagrangian 
with an one-particle term and an effective electron-spin coupling of form 



effective spin operator, the leading order coupling terms between spins and density /current 
fluctuations can be derived and are represented in the Feynman Diagram shown in Fig. 
(la), where the solid lines are electron propagators. In real space-time, the diagrams are 




(12) 



where ip = (c-f,C|) is a 2-component spinor where c a 's are electron operators and S is an 
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represented by expressions of form 

r^(x,x'_,x' + ; S) = ^ G a (x - x'_)j^(x)G a (x' + - x)S v (x'J)r^ a ,G a '(5x)S ,/ '(x' + )T^ a (13) 

aa'vu' 

= J2 G *( X - x'S]»(x)G a (x' + - x) [G- a (5x)(S x (x'_)S x (x' + ) + S y (x'_)S y (x' + )) 

IT 

+ G a (5x)S z (x'_)S z (x' + ) - i(a)G^(5x) (S x (x'_)S y (x' + ) - S y (x'_)S x (x' + ))] 

where j^fi = 0, 1,2) is the electron current operator and x = (x,t), = x' — (+)8x/2. 

Assuming that the electrons have a gapped spectrum (insulator), the corresponding Green's 
function G a (x) is short-ranged and the contributions mainly come from small 8x region. 
Therefore we can expand G a (x' ± 5x/2 — x) ~ G a (x' — x) ± (5x/2)d x G a (x' — x) + .., 
S u (x' ± Sx/2) ~ S u (x') ± (8x/2).d xl S v + .., etc. in Eq. to derive the leading order 
spin-density(current) coupling terms in the insulating state in the continuum limit. A cor- 
responding expansion for metallic ferromagnetic states has been developed previously^. By 
keeping two sites per unit cell, this procedure can be extended to derive the correction 
to optical conductivity in the antiferromagnetically ordered state in the Hubbard model, 
which is a competing state to the spin-liquid state observed in the organic compound 
k - (BEDT - TTF) 2 Cu 2 (CN) 3 10 . In this case, G(x) -> G ab (x) and S(x) -> S a (x) = 
m(x) + (— l) a n(x), where a, b = A,B are sublattice indices, in and n represent magneti- 
zation and staggered magnetization fluctuations, respectively. The low energy contribution 
to optical conductivity is dominated by coupling of density fluctuations to two spinwave 
process represented by coupling to n fields. After some algebra, we obtain in the small 
wave- vector limit, 

T (q, o>; q', Q;q — q',ui — Q; S) ~ u(q.(f)n((f , fl).n(q — <f,uj — Q). (14a) 

We have assumed that the antiferromagnetic state is described by usual mean-field theory 
with non-zero staggered magnetization < m >. The corresponding correction to proper 
density density response function (Fig. (lb)) is given by 

, ,„ , 1 V- |r (0.^;,,'.iW.^-fi;S)| 2 

SM ~w^ a w^m^w^?' (14b) 

where c m ~ U < m > is the spinwave velocity derived from the mean-field theory. Evaluating 
the integral, we find that the correction to optical conductivity is 

Safa) ~ e 2 ( — ) d+2 , (14c) 

Cm 
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for 00 << U < m >, where d is the dimension. We have assumed t <~ U in our calculation. 
Notice that the optically conductivity is enhanced in the spin-liquid state compared with the 
antiferromagnetically ordered state, in agreement with what is observed experimentally^. 

In conclusion, we have shown that gapless spinons in a spin liquid state gives rise to a 
power-law optical absorption inside the Mott gap which is larger than that expected for two 
spin wave absorption in a Neel ordered insulator. Recent experiment has reported the sur- 
prising finding that the low temperature optical absorption in k,—{BEDT—TTF)2Cu2{CN)^ 
is larger than another compound k — (BEDT — TT F) 2 Cu[N (C N) 2 ]Cl which exhibit Neel 
ordering but is " closer" to the Mott transition in that it has a smaller Mott gapi^. Our result 
gives a natural explanation of this puzzle. We believe that power-law absoption, especially 
if it can be observed in a large gap insulator such as the Kagome system, is strong evidence 
for the existence of gapless spinons and gauge fields. 
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